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Abstract 

The analytical perturbative approach, if taken to the limit of its applicability, 
allows one to predict an energy independent limit for the one-particle invariant 
density in QCD jets E$- at very small momenta p. This is a direct conse- 
quence of the colour coherence in soft gluon branching. The existing data 
on the charged and identified particle inclusive spectra follow this prediction 
surprisingly well. Further tests of the perturbatively based picture in the soft 
region are discussed. 
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1. Introduction 



One of the most striking predictions of the perturbative approach to QCD jet physics |j]J 
is the depletion of the soft particle production and the resulting approximately Gaussian 
shape of the inclusive distribution in the variable £ = log Ej et /E for particles with energy 
E in a jet of energy Ej et (the so-called "hump-backed plateau" ) |||4]] . Due to the coherence 
of the gluon radiation it is not the softest partons but those with intermediate energies 
(E ~ Ej£~ 0A ) which multiply most effectively in QCD cascades. We recall that the 
gluons of long wave length are emitted by the total colour current which is independent 
of the internal structure of the jet and is conserved when the partons split. Applying 
the hypothesis of the Local Parton Hadron duality (LPHD)||, one then comes to the 
conclusion that the hadron spectrum at low momentum p should be nearly independent 
of the jet energy E jet ^. 

It is a very important question to understand whether there is a smooth transition 
between the purely perturbative regime and the soft momentum region. It looks quite 
intriguing that the shape of the measured charged particle energy spectra in e + e~ anni- 
hilation appears to be surprisingly close over the whole momentum range (down to small 
momenta of a few hundred MeV) to the perturbative predictions based on the Modified 
Leading Log Approximation (MLLA) ||||[7||| . Moreover, the momentum spectra in the 
variable logp, measured by the TASSO and OPAL Collaborations [|Kj provided first 



evidence in favour of the scaling behaviour at small momenta: as discussed in |l Ij these 
data could be considered as a confirmation of the basic ideas of QCD coherence and 
LPHD. A further confirmation of this picture has been provided recently in the study of 
the invariant particle density Edn/d 3 p at low momenta and also in the low cms energy 
region. The spectra are found to be in a good agreement with the scaling behaviour and 
with the analytical perturbative expectations which become sensitive to the running of 
the coupling a s at small scales [|1^ . 



In this paper we focus on the soft limit of the parton spectra, and how its scaling 
behaviour arises in the analytical calculations for the coherent parton cascade. We con- 
sider the invaraint density Edn/d 3 p = dn/dyd 2 pT in the limit of vanishing rapidity y and 
transverse momentum px or, equivalently, for vanishing momentum |p| = p, i.e., 

Io= Hm 4 = ilim^ (1) 

y^o,p T ^o a 6 p 2 p^o d A p 

where the factor 1/2 takes into account that both hemisphere are added in the limit p — > 
0. If the dual description of hadronic and partonic final states is really adequate down to 
very small momenta, the finite, energy independent limit of the invariant hadronic density 
Iq is expected. 

This will be examined for both charged and identified hadrons in the full energy range 
explored so far in e + e~ annihilation. A possible growth of the spectrum in the soft region 
could indicate that either coherence or the local duality (or both) breaks down. Since 
colour coherence is a general property of QCD as a gauge theory, it is the LPHD concept 
that is tested in the measurements of the soft hadron distributions. Therefore, such 
studies could provide one with additional information on the physics of confinement. 

The scaling properties of the invariant cross section or density have been discussed in 
the early seventies in the context of the parton model [|13[ and Feynman-scaling |14| (see 
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e.g. Note that scaling in the parton model refers to the pr-integrated spectra. In the 

QCD bremsstrahlung picture of multiparticle production discussed here the p^-integrated 
central rapidity plateau rises with cms energy, and scaling behaviour is retained only in 
the large wavelength limit ([!]). 

It is important to test further the QCD origin of this scaling behaviour. This can 
be accomplished by studying the soft radiation from the emitters ("parton antennae") 
corresponding to different effective colour charges. The good scaling property of the 
density Iq in e + e~ annihilation suggests taking this number as a standard scale in the 
comparison with other processes. Some of the results of this paper and further details 
have also been presented elsewhere |T7| . 



2. Particle production in the soft limit in DLA and MLLA 

DLA predictions. 

We consider first the analytical predictions for the energy spectrum of partons near the soft 
limit. The asymptotic behaviour of the energy spectrum is obtained from the Double Log 
Approximation (DLA) in which energy conservation is neglected and only the leading 
singularities in the parton splitting functions are kept. The evolution equation of the 
inclusive energy distribution of partons p originating from a primary parton A is given 
by§: 

Y) = 5 A 5(£) + f d? n dyi^ (y>)D^>, y>) . (2) 
Jo Jo JMc 

Here we have used the logarithmic variables £ = log(l/a;) = \og(Q/E) and Y = \og(Q/Qo) 
with E the particle energy and Q the jet virtuality (Q = PQ for a jet of primary mo- 
mentum P and half opening angle 0); Ca is the respective colour factor, i.e., C g = N c 
and C q = Cp; 7 denotes the anomalous dimension of multiplicity and is related to the 
QCD running coupling by = 4N c a s /27t or 7q = (3 2 / log(p T /A) with /3 2 = 4N c /b, 
b = (11N C — 2nj)/3; A is the QCD-scale and Nc and n/ are the number of colours and 
of flavours respectively. The shower evolution is cut off by Qo, such that the transverse 
momentum pp > Qo- 

As can be seen from (Q) the ^-spectrum vanishes in the soft limit £ — ► Y (E — > Qo). 
The behaviour near this limit can be found by the iterative solution of (|2|). With two 
iterations one obtains 
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/? 2 / y -^rlog(l + ^)log(l + 4 x ) 
log(l + ^) 



+. 



(3) 

where A = logQo/A. The second term of order (3 2 corresponds to the emission of a single 
gluon and yields the leading contribution for E — > Q . It is proportional to the colour 
charge factor of the primary parton. Furthermore, this term does not depend on the cms 
energy, contrary to the higher order terms which provide the rise of the spectrum for large 
E with increasing yfs. The spectrum vanishes at E — > Qo as 

D 9 A {^ Y) ~ Y - e ~ log E/Qo ~E-Qo. (4) 

An exact solution is known for fixed a s [^] with similar behaviour near the boundaries as 
above in (j^) and (J|). 
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MLLA predictions. 

A considerable improvement in accuracy for nonasymptotic energies is obtained in the 
application of the MLLA |6|,0,0 which takes into account systematically all the corrections 
of relative order ^fcT s based on the evolution equation which includes energy conservation 
and the exact form of the parton splitting functions. The leading high energy MLLA 
contribution obeys a differential equation (see also [@,0, with D(£, Y, A) = D 9 g ) 

+ JU* loD(^Y,X) = -a — + — — D(C,Y,X)\ (5) 



d£dY dY 2 ,u w ' ' \di dY)\ 2tt 

where a = y^c + The limit a = corresponds to the DLA as in (0). In case of fixed 
a,* the exact MLLA solution can be obtained jL^] (see also ||) by multiplying the DLA 



result with the factor exp(— a^^Y — £)/4iVc). Similarly, for running a s one can make an 
ansatz 



D(Z,Y,\) = D^,Y\)\DLAexp 



-a f \o(y)/(4N c )dy 



(6) 



which solves Eq. ([5]) approximately in the region around £ ~ Y.f\ Note that the exponent 

in (||) can be written as — | log^^y^ and behaves near the boundary like — | , so the 

slope of the distribution varies with energy like the coupling a s ~ 7q = 4Nc/[b(Y + A)] 
contrary to the case with fixed coupling. 

Phase-space effects in the soft region. 

While the limiting behaviour of the partonic energy spectrum in the soft region follows 
from the general principle of colour coherence the detailed form of the observable hadronic 
spectrum is predicted uniquely from the LPHD hypothesis only for E ~ p ^> Q , but 
not near the kinematical boundary because of the sensitivity of the spectrum to the cut- 
off procedure. In the MLLA the partons are treated as massless with energy E = p > 
Pt > Qo, so £ < Y. Experimental hadronic spectra are usually presented as function of 
momenta p or £ p = log(l/x p ) which is not limited from above. The same kinematic limit 
for partons and hadrons is obtained if the hadronic mass rrih and the partonic pt cut-off 
Qo are taken the same. 

For the relation between parton and hadron distributions one may require that the 
invariant density Edn/d 3 p of hadrons approaches a constant limit for p — > as is observed 
experimentally. For the spectra which vanish linearly as in (f|) this is achieved by relating 
the hadron and parton spectra in a single A-jet as @,|TJ] 

Eh — -j = K h E p — = K h D A (£ E , y ) {<) 

dp h dp p 



with E h = Jp\ + Qo = E P = E > Q and ^ E = C, = logQ/E, where K h is a normalization 
parameter. If hadrons from both hemispheres are added, Kh should be replaced by 2K^. 



4 When eq. (Q) is inserted into eq. (g) all terms cancel except for the small difference 
(dD DLA /dY) exp[— a ^(y) / {AN c )dy\[ aa ^ — aa ^' ) ] 1 which vanishes at £ = Y. Exponential damping 
factors of this type are known to take into account single log corrections [Oy. 
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Then, indeed, for hadrons the invariant density Edn/d 3 p = K h D A (C i ,Y)/4:n(El — Qq) 
approaches the finite limit as in ([!]) 

T K c A p 



8nN c \Q 2 



In the fixed a s limit [3 2 /\ is replaced by 70 and Iq ~ 1 / Qo- With prescription (0) the 
moments of the full energy spectrum D(£,Y) are well described by the MLLA formulae 



and with Q = 270 MeV in a wide energy range (12[ . 

The relation (|7|) is not unique, however. We found that the alternative prescription 
based on phase space arguments, dn/d^ p ~ (p/ E) 3 D(!;e,Y) (see e.g. [ITglBDJ): 



E *wr K ^) D3 ^' Y) (9) 

works well in e + e~ annihilation for charged pions (Kh = K n ) in the whole energy region 
and for charged particles (Kh = K c h) in the low energy region. In this application, 
including both hemispheres, D 9 A in eq. @ is replaced by 2-4/9-D^, using for D 9 g the MLLA 
limiting spectrum PJT8[] with Qq = A = 138 MeV. The low cut-off mass is plausible in 
this region, which is dominated by pions. 
Behaviour of analytical results 

In order to illustrate the above analytical results we compare in Fig. 1 the various approx- 
imations (in case of running a s ) for low particle energies using the relation (p]) between 
parton and hadron spectra: the DLA result in the approximation @ with and without 
MLLA exponential factor (||) for A = 0.01 and also the "limiting spectrum", which solves 
the MLLA equations for A = in the full kinematic region, except very close to the 
boundary £ = Y where it stays finite. The normalization of the limiting spectrum is as 
in fits to e + e~ annihilation; the normalizations of the DLA and MLLA curves are chosen 
to approach the limiting spectrum for energies E > 0.5 GeV. 

These results show that in all cases considered, DLA and MLLA for both fixed and 
running a s , the single particle invariant density approaches an energy independent value in 
the soft limit £ — > Y (E — > Qq). This originates from the soft gluon emission contribution 
of order a s which is determined by the total colour charge of the primary partons due 
to the colour coherence. In this limit the MLLA converges towards the DLA as the 
energy conservation constraints are unimportant and the parton splitting functions are 
only probed for very small fractional momenta. 



3. Discussion of experimental results 



Inclusive charged particle spectra 

Fig. 2 shows the charged particle invariant density in e + e~ annihilation, Edn/d 3 p, as a 
function of the particle energy E at different cms energies ranging from 3 GeV up to 
LEP-1.5 cms energy (133 Ge V) g|TD|,g2],g3] . The effective mass value of Q = 270 MeV 



has been used in the kinematical relations. As shown in [0 this allows a reasonably good 
description of the moments of the energy spectra by the MLLA limiting formulae over a 
large energy interval. A fit using eq. (|9|) is shown elsewhere (see last reference in |JL7|). 
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It is remarkable that the data from all cms energies tend to converge in the soft limit: 
we find 2J ~ 6-8 GeV" 2 (using Q = 270 MeV) and 2J ~ 4-6 GeV" 2 (using Q = 138 
MeV). Inspecting the soft limit more closely, the LEP data seem to tend to a limiting 
value larger by about 20% as compared to the lower energy data. This may be well due to 
the overall systematic effect in the relative normalization of the different experiments. In 
this respect, it is instructive to recall that both TASSO and the LEP Collaborations have 
collected their data at different cms energies with the same accelerators and detectors, 
thus avoiding within their data samples any problems due to the relative normalization; 
in these samples the scaling works well. 

A possible source of scaling violations are weak decays (for example kaons, heavy 
quarks). Particles from such decays (from the new "antenna") would add incoherently 
to the particles produced from the primary quarks and thereby could yield a rise of the 
soft particle spectrum with increasing energy. It would be interesting to disentangle these 
possibilities experimentally and to find out the size of primary scaling violation, if any. 

The data also compare well with the MLLA analytical predictions, not only concerning 
the energy independence of the soft limit Jo but also the considerable energy dependence of 
the initial slope which is largely due to the running of a s in the MLLA damping factor (|6|). 



Identified particles ' spectra 

Let us consider now the invariant density for tt, K and p. In this case, the simplest 
assumption is to identify with the mass scale the particle mass itself. Fig. 3 shows the 
invariant cross section Edn/d 3 p as a function of the particle energy E for charged pions, 
charged kaons and protons, as derived from the inclusive momentum spectra measured 



at cms energies from 1.6 GeV to 91 GeV[ p3|j24]1 . In all cases the data tend towards a 
common limit for E — > m^ (p — > 0). 

Note that the scale for I in eq. (H) is given not by A~ 2 , but by Qq 2 , which can be 
related to an effective particle mass. Indeed, Iq shows a clear mass dependence (2/q 
= 5-8 GeV" 2 , 21* = 0.3-0.4 GeV~ 2 , 21$ = 0.08-0.10 GeV" 2 ). This is consistent with 
Iq ~ m^ 2 (after normalization to Iq , one would predict Iq = 0.37-0.61 GeV -2 , Iq = 
0.10-0.17 GeV -2 ). Inclusion of the logarithmic correction 1/A in eq. (§) would, however, 
require rather small values of A ~ 100 MeV. The measured values of Iq are also in good 
agreement with the expectations following from eq. @ and the finite limit at E — > Qq 
of the limiting spectrum. Alternatively, one can use in the fits effective scales Qq which 
deviate from the respective hadron masses; this is also suggested by the fits to the peak 
region of the £-distributionF]. This intriguing subject certainly deserves further detailed 
studies. 



4. Further tests of the dual description of parton and hadron cascades in 
e + e annihilation. 

It is suggestive to relate the observed behaviour of particle spectra for vanishing momenta 
to the expected behaviour of soft gluon emission off the primary partons. The particle rate 
Jo in this limit (see eq. @) cannot be directly predicted, as it depends on the normalization 

5 For example, the maximum £* of the distribution does not follow the expectations (£* ~ 
l/21og£ jet /Q in DLA) with Q = m h ^. 
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factor and the cut-off parameter Q . However, the remarkable scaling behaviour of 1$ 
in e + e~ annihilation provides one with a standard scale for the comparison with other 
processes. It would be very interesting to study whether the soft particle production is 
universal, i.e., a purely hadronic quantity or whether the intensity I indeed depends on 
the colour topology of the primary active partons in the collisions process. This could 
provide one with a direct support of the dual description of soft particle production in 
terms of the QCD bremsstrahlung. 

In what follows we consider some possibilities to test this hypothesis further and dis- 
cuss what one may expect from the comparison of the e + e~ reaction with other collision 
processes. 



gg final state 

A direct test of the relevance of QCD coherence for the soft production limit is the 
comparison of the qq with the gg colour singlet final state. In the gg final state the 
same line of argument applies as above for the qq final state, but the colour charge 
factor is increased by N c /Cp = 9/4 in eq. (§) implying roughly a doubling of the soft 
radiation intensity Jo- This increased radiation has been suggested originally for the global 
multiplicity in gluon jets[26[; in the real life experiments various effects keep the observed 
increase considerably below this limit (see e.g. ||). As the soft radiation is not strongly 
influenced by energy-momentum constraints (see the DLA-MLLA comparison above) the 
effects of the different colour charges could be more pronounced in this case. 

An approximate realization of a colour octet antenna is possible in e + e~ — > qqg with 
the gluon recoiling against a quasi-collinear qq pair[^7| (for a recent study of such events, 
|28|). The production rate of soft particles in such events is expected to be increased 



sec 



if it really is related to the effective colour charge of the primary partons. 



3- jet events 

To be more quantitative, we consider the soft radiation in 3-jet qqg events of arbitrary 
jet orientation into a cone perpendicular to the production plane and compare it to the 
radiation into the same cone in a 2-jet qq event again perpendicular to the primary qq 
directions. The analysis of these configurations avoids the integration over the kx > Qo 
boundaries along the jets. 

We restrict ourselves to calculations of the soft gluon bremsstrahlung of order a s . 
From the experience above with two-jet events this contribution dominates in the soft limit 
whereas higher order contributions take over with increasing particle energy. Calculations 
of this type have been successfully applied to the string/drag phenomena [ p9] , |30H which refer 



to the particle multiplicity flows projected onto the 3-jet production plane. 

Let us consider first the soft radiation into arbitrary direction n from a qq antenna 
pointing in directions n, and Hj [30 : 



dN ^d^J^WM^) W^{n) = 2C F {ij) (10) 

with (ij) = dij/faidj), ciij = (1 — Hifij) and a; = (1 — nn^). Such an antenna is realized, 
for example, in a qcpi event, and can be obtained by the appropriate Lorentz boost from 
the qq rest frame. The soft gluon radiation in a qqg event is given as in ( |I~0"D but with the 
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angular factor 

W™{n) = N C [(U) + (T=) - ^CR] (11) 

where (+, — , 1) refer to (q, q, g). 

For the radiation perpendicular to the primary partons (ij ) = = 1 — cos 0jj , with 
relative angles 0^- between the primary partons % and j. For 2-jet events of arbitrary 
orientation one obtains in this case 

Wf(6+_) = 2CV(1 - cos 0+_) . (12) 

Correspondingly the ratio of the soft particle yield in 3-jet events to that of 2-jet 
events in their own rest frame (Wj q (ir) = 4CV) is given by 

i2± = "J^T = - cos @1 + - cos 01 - " iv^ (1 " C0S @+ - )] (13) 

J_ F c 

It is also interesting to note the difference of this prediction to the large Nc approximation 
in which the qqg event is treated as a superposition of two qq dipoles (see, e.g., f29| ). In 



this case the last term in eq. fll3D drops out and Cp = (Nq — 1)/2Nc — Nc/2. This yields 

R ± ~ - [2 - cos 0i+ - cos 0!_] (large A^ c ) (14) 

Predictions from these formulae for R± are presented in Table (|]) for various relative 
angles 0^. Note, in particular, the limiting cases R± = 1 for soft or collinear primary 

gluon emission and the proper gg limit for the parallel qq (0_| = 0) configuration, as 

expected. 

The role of the large- N c limit can be investigated also by studying the production 
rate in 3-jet events normalized to the sum of rates from the corresponding 2-jet events 
(dipoles) with opening angle 0i+ and ©i_ respectively: 

x = divf(© 1+ ) + divf(0 1 _) (15) 



Nl , 1 l-cos(0i+ + 



Nl - 1 V Nq 2 — cos 0i+ - cos 0i_ 

This ratio measures directly the deviation from the large-TVc limit R±= 1 and thereby 
from the gg-dipole approximation. This approximation is not necessarily limited towards 
soft particle production. For the simple case of Mercedes- like events (0i+ = 0i_ = ©+-) 
one obtains R± = 17/16 = 1.06. The 2-jet rates for relative angle 0^ which appear in 
the denominator of eq. (|1^) could be found experimentally from the corresponding q<pi 
final states. 

Let us list a few further results: 

a) A particularly simple situation is met for Mercedes- type events where no jet identifi- 
cation is necessary for the above measurements. 

b) The large angle radiation is independent of the mass of the quark (for ^> mQ/Ej et ). 

c) The above predictions for the ratios R±, R± are derived for the soft particles according 



to the bremsstrahlung formula (]10). One may also consider the particle flow integrated 
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Ri 



R ± (large N ( 



c 



71 



e 1+ = vr - 9x_ 

(collinear or soft gluons) 

e 1+ = e+_ 
e 1+ = ©+_ 
(->.. (-), - .,. 

(Mercedes) 

e 1+ = 6i_ = vr 

(qq antiparallel to g) 



7T 



1.21 
1.42 
1.59 



^ = 2.25 

Of 



1.18 
1.35 
1.5 



Table 1: Prediction for the ratio R± = dNf 9 /dNf from fll3|) and its large- Nc- 
approximation (|j~4|) for different configurations of the qqg events (Oi+ = Q gq , ©i_ = Q g g, 
0+_ = 2tt- 6i+ - 6i_). 



over momentum, as in the discussion of the string/drag effect. In this case one has to 
include all higher order contributions which take into account the fact that the soft gluon 
is part of a jet generated from a primary parton. Then the angular flow dN/dQ^ is given 
by the product of the radiation factor W{n) and a "cascading factor" |3(J. For the ratio 
of multiplicity flows one obtains the same predictions , fllBD as the cascading factors 
cancel. It will be interesting to see to what extent the predicted angular dependence 
for both quantities - the multiplicity flows and the soft particle yields - are satisfied 
experimentally. 

d) The similarity of particle flows in 3-jet and 2-jet events with corresponding angles 0i+, 
0i_ , as expected in the large- Nc approximation, should also apply to further details of 
the final state such as the particle ratios. Since the Lorentz transformation along the 
boost direction produces a larger drag for heavier particles, one may expect that after the 
boost the K/tc and p/n ratios for soft particles decrease and then the same is true for soft 
particles in 3-jet events in comparison to 2-jet events in their rest frame^. 



5. Extension to different reactions 

We consider here processes with primary hadrons or photons with dominantly 2-jet final 
states. This includes semihard processes which are initiated by partonic 2-body scatter- 
ings. In case of quark exchange the two outgoing jets originate from colour triplet charges 
and Jo should be as in e + e~ annihilation; in case of gluon exchange Iq should be about 
twice as large (Nc/Cp) as in a gg jet system. In 2-jet events the limiting density I may 

6 This expectation has been verified in the JETSET Monte Carlo[|3l|]. We thank T. Sjostrand for 
providing us with this information. 
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be most conveniently determined as limit of dn/dycPpT at y ~ for p? — » o. Final states 
with several well separated jets can be treated in analogy to the e + e~ — > 3-jet events 
discussed before. 

fa) owarJ; exchange processes 

A simple example is the process 77 — > gg — > 2-jets which dominates if either the virtu- 
ality Q 2 of the initial photon or the scattering angle photon-jet is sufficiently large (see, 
e.g. [0]). In this case the soft production intensity Jo should be the same as in e + e~ 
annihilation. 

Another example is deep inelastic scattering at large Q 2 . In this case the current 
fragments in the Breit frame are expected to have the same characteristics as the quark 
fragments in one hemisphere of e + e~ annihilation (for a QCD analysis, see ||33|| ), and this 
is indeed observed for not too small Q 2 133],|35] . Therefore one expects again the same limit 



Jo as in e + e annihilation. 

(b) gluon exchange processes 

These are expected to dominate for virtual 7y7 or jvP scattering (Q 2 > few GeV 2 ) 
through the photon gluon fusion subprocess 7^ — > qq at small Bjorken x. In this case (for 
a fast qq pair in one hemisphere with small opening angle) the particle jets emerging in 
opposite directions are generated by an effective colour octet emitter and the soft intensity 
Jo should be about twice as large as in e + e~ annihilation. 

Another example is hadron-hadron collision with a particle or jet at moderate pt (> 
1-2 GeV) at small angles so that the overall 2-jet structure is maintained. This type 
of scattering is dominated by gluon exchange and the same increase of Iq is therefore 
expected. 

(c) soft collisions (minimum bias events) 

These processes (with initial hadrons or real photons) are not so well understood theo- 
retically as the hard ones but it might be plausible to extrapolate the gluon exchange 



process towards small pr[3o]. Experimental data at ISR energies, however, do not fol- 
low the expectation of a doubling of Jo in comparison to e + e~ annihilation, rather the 
intensities are found to be similar [57|]. On the other hand, Jq roughly doubles when going 



from v/i ~ 20 GeV @,|3§ to ^/s = 900 GeV at the collider Q. If additional incoher- 
ent sources like weak decays can be excluded, such behaviour could indicate the growing 
importance of one-gluon exchange expected from the perturbative picture. Then a satu- 
ration at Ig h j Iq +e ~ 2 is expected and no additional increase for a semihard pt trigger. 
A rise of Iq 11 could also result from incoherent multiple collisions of partons (e.g. f41| ) 
which recently has also been postulated for jp collisions ||42||. 



It should be noted that events from all three classes a, b and c can occur in the same 
reaction depending on the external kinematical constraints (triggers). For example, if 
Q 2 for fixed hadronic energy W in jp or 77 collisions is decreased from large Q 2 ~ W 2 
towards zero one changes from (a) to (b) and finally to (c). Also, a dependence on the 
reference frame is expected, case (a) in Breit frame, case (b) typically in cms frame for 
small Bjorken x. 

Rapidity dependence of Iq. 

The frame dependence can be studied conveniently in these exchange processes by con- 
sidering Iq (y) = dy j2 pT \ PT ^ as a function of rapidity y, measured, say, in the cms frame. 
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The Breit frame is reached by a Lorentz transformation in the direction of the incom- 
ing photon, and corresponds to a different rapidity, y Breit, m the cms. According to the 
above discussion, one would then expect for the processes with virtual photons at large 
Q 2 in general a "quark plateau" of Io(y) in the current region of length Ay pa \og(Q/m) 
near ysreit and a transition to a "gluon plateau" in the complementary region of length 
Ay pa \og(W 2 /Qm) where m is a typical particle mass. The existence of a plateau cor- 
responds to the energy independence of I as seen in e + e~ annihilation. So, if gluon 
exchange occurs with sufficient hardness in the process, Iq should develop a step like be- 
haviour but never become larger than twice the Jo value in e + e~ annihilation. The height 
of the "plateau" is a direct indicator of the underlying exchange process (quark or gluon 
like). 

6. Conclusions 

The analytical perturbative approach to multiparticle production, based on the Modified 
Leading Log Approximation and Local Parton Hadron Duality (LPHD) has proven to be 
successful in the description of various inclusive characteristics of jets. It is of importance 
to investigate the limitations of this picture, in particular in the soft region where non- 
perturbative effects are expected to occur. 

It is remarkable that the intensity of the soft hadron production follows to a good 
approximation a scaling law in a range of two orders of magnitude in the cms energy of 
e + e~ annihilation (1.6-140 GeV). Such a scaling law is derived analytically for the soft 
gluons in the jet and follows directly from the coherence of the soft gluon radiation from 
all emitters. It appears that the production of hadrons which is known to proceed through 
many resonance channels nevertheless can be simulated in the average through a parton 
cascade down to a small scale of a few 100 MeV as suggested by LPHD. The scaling 
behaviour down to the very low cms energy of 1.6 GeV can only be explained within a 
parton model description if the cut-off scale Q is well below 1 GeV. 

It will be interesting to investigate the suggested scaling property (1) and its possible 
violation further in the same experiment to avoid systematic effects; this seems to be 
feasible at LEP and HERA for the simplest processes. In particular, the effect of weak 
decays on scaling violations should be clarified. The sensitivity of the soft particle pro- 
duction to the effective colour charge of the primary emitters can be tested through the 
transverse production rates in multijet events. The soft radiation in e + e~ annihilation 
can be used as a standard scale in the comparison of various processes. In this way the 
soft particle production can be a sensitive probe of the underlying partonic process and 
the contributions from incoherent sources. 
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■ MLLA (A = 0.01) 

- DLA (A = 0.01) 

Limiting Spectrum (A = 0) 
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Figure 1: Invariant density Edn/d 3 p as a function of the particle energy E for Qq = 
270 MeV. Predictions at cms energies of y/s = 3 GeV (lower three curves) and 91 GeV 
(upper three curves) using Edn/d 3 p = 2 ■ A/QK h D 9 g (£e)/[47t(.E 2 — Ql)] with 7J| computed 
in MLLA (eqs. (§,0), DLA (eq. (||)) and the Limiting Spectrum (normalization K h = 0.45, 
0.45 and 1.125 respectively). 
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E[GeV] 



Figure 2: Invariant density Edn/d 3 p of charged particles in e + e annihilation as a function 
of the particle energy E = \jp 2 + Ql at Q = 270 MeV. Data points at various cms 
energies from SLAC, TASSO and TOPAZ Collaborations, LEP-1 and LEP-1.5f23)||[n| 
are compared to MLLA predictions (normalization as in Fig. 1 with A = 0.01, = 0.45). 
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Figure 3: a) Invariant density Edn/ d 3 p of charged pions in e + e annihilation as a function 
of the particle energy E at Qo = 138 MeV at various cms energies; b) the same as in a), 
but for charged kaons (with Q = 494 MeV) and protons (with Q = 938 MeV) fUS 
The kaon distribution is multiplied by a factor 10 for the sake of clarity. 
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